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ABSTRACT:  Composite  materials  have  complex  failure  modes  that 
include  delamination,  fiber  debonding  and  breakage,  and  matrix  mi¬ 
crocracking,  In  this  paper,  the  influence  of  these  damage  modes  on 
the  failure  of  the  short-beam  three-point  bend  test  is  investigated  for 
a  composite  with  a  quasi-isotropic  layup.  Failure  is  found  to  initiate 
in  a  region  near  the  point  of  application  of  the  load,  a  location  where 
classical-type  analytical  descriptions  of  specimen  behavior  are  unre¬ 
liable.  Furthermore,  the  locations  of  failure  show  little  reproducibil¬ 
ity.  Observed  fracture  behavior  is  explained  in  terms  ot  the  overall 
stress  state  in  the  beam  before  fracture,  and  failure  is  predicted  from 
the  stress  map.  using  the  maximum  strain  criterion. 

KEY  WORDS;  composites,  delamination,  fractography.  finite  ele¬ 
ment  analysis,  failure  analysis 

The  three-point  short-beam  bending  test  (ASTM  Test  for  Ap¬ 
parent  Interlaminar  Shear  Strength  of  Parallel  Fiber  Composites 
by  Short-Beam  Method  D  2344),  as  shown  in  Fig.  1,  has  long 
been  used  to  quantify  the  shear  strength  of  laminated  composite 
materials.  The  interlaminar  shear  stress  at  failure  is  then  cal¬ 
culated  using  the  equation 


where 

T,  =  interlaminar  shear  stress, 

P  =  applied  load  at  failure,  and 
A  =  cross-sectional  area  of  the  specimen. 

The  standard  procedure  for  the  short-beam  shear  test  is  in¬ 
tended  for  use  with  unidirectional  composites  with  fibers  oriented 
along  the  length  of  the  specimen.  Equation  1  is  based  on  a 
parabolic  through-the-thickness  shear  stress  distribution  with  a 
maximum  at  the  beam  midplane.  This  simple  strength  of  ma¬ 
terials  analysis  neglects  the  effects  of  stress  concentrations  and 
material  anisotropy. 
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This  test  is  frequently  used  to  screen  materials  based  on  their 
interlaminar  shear  strength  However,  there  is  rof-  '^x-n  .jse 
of  the  test  relies  upon  incorrect  assumptions  concerning  the  na¬ 
ture  of  the  failure  process. 

Whitney  and  Browning,  in  Ref  1.  in  examining  the  failure  of 
unidirectional  short-beam  specimens,  discovered  that  pure  Mode 
II  failure  along  the  specimen  midplane  was  not  the  mode  of 
fracture.  Rather,  compressive  stresses  in  the  regions  of  high  in¬ 
terlaminar  shear  stress  were  found  to  suppress  interlaminar-shear 
failure  modes,  and  initial  damage  as  a  result  of  vertical  cracking 
was  found  to  precede  final  failure  caused  by  large-scale  delam¬ 
ination.  If  the  failure  process  is  this  complex,  then  the  applied 
load  at  failure  does  not  depend  solely  on  the  interlaminar  shear 
strength,  making  the  test  useless  for  its  intended  purpose. 

The  elastoplastic  finite  element  analysis  of  Berg  et  al.  [2]  showed 
extremely  high  shear  stresses  near  the  load  nose,  and  the  finite 
difference  analysis  of  Sandorff  (i)  showed  that  St.  Venanfs  prin¬ 
ciple  is  not  satisfied  in  an  orthotropic  beam  of  low  span-to-depth 
ratio,  with  the  result  that  the  stress  concentration  effects  near 
the  load  nose  and  the  supports  are  never  dissipated. 

Because  of  this  uncertainty,  there  is  a  need  for  reexamination 
of  the  short-beam  three-point  bend  specimen  for  specific  com¬ 
posite  configurations.  Such  an  examination  must  determine 
whether  a  simple  relationship  between  loading  and  fracture  ex¬ 
ists,  and  if  not.  what  fracture  processes  are  taking  place.  The 
current  effort  answers  these  questions  for  a  graphite/epoxy  lam¬ 
inate  with  a  quasi-isotropic  layup. 


FIG.  I  —  The  ihrce  poini  shori-heam  lesi. 
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General  Failure  Analysis 

Specimen  Design  and  Testing 

The  subject  of  the  investigation  was  a  short-beam  specimen 
with  dimensions  given  in  Table  1.  A  sample  panel  was  fabricated 
and  inspected  for  quality,  using  standard  aerospace  industry  prac¬ 
tice,  befoie  sectioning  into  specimens.  Two  fiberglass  plies,  each 
of  thickness  0.457  mm  (0.018  in.),  were  bonded  to  the  upper 
and  lower  surfaces  to  minimize  the  bearing  effect  at  the  supports 
and  the  applied  load  location. 

The  material  properties  for  the  graphite/epoxy  and  the  fiber¬ 
glass  plies  are  given  in  Table  2.  The  material  was  an  AS4/3502 
graphite/epoxy  system,  and  the  layup  was  [+45/0/90]„,.  Failure 
load  data  for  four  specimens  tested  to  failure  in  three-point  bend¬ 
ing  are  given  in  Table  3.  The  shear  strength  predictions  in  the 
table  are  based  on  Eq  1.  The  failed  specimens  were  sectioned 
with  a  silicon  carbide  saw,  and  the  sections  were  analyzed  using 
a  scanning  electron  microscope  (SEM). 

Numerical  Methods 

Numerical  analysis  of  the  specimen  concentrated  on  the  phe¬ 
nomenon  of  delamination.  Delamination  analysis  can  be  based 
on  two  approaches;  the  strain  energy  release  rate  approach  and 
the  interlaminar  stress  approach.  The  interlaminar  stresses  de¬ 
velop  at  interfaces  between  plies  to  keep  the  laminate  in  a  state 
of  equilibrium  and  kinematic  compatibility.  Delamination  occurs 
when  these  stresses  reach  the  interlaminar  strength  of  the  ma¬ 
terial.  The  strain  energy  release  rate  approach  is  based  on  the 
actual  process  of  fracture,  rather  than  the  strength  concept.  De¬ 
lamination  can  propagate  when  the  strain  energy  release  rate  at 


TABLE  1 — Dimensions  of  quasi- 
isotropic  short-heam  specimen,  mm 
(in.)  (.see  Fig.  I  for  nomenclature). 


L 

24.4(1  (l.(KI) 

c 

15.S75  (0.62.S) 

t,..h 

l).4.S7  (O.OIX) 

7.188  (0.283) 

B 

9.,S2.S  (0.375) 

TABLE  2  —  Properties  of  specimen  constituent  materials. 


Constant 

AS4:'3502 

Fiberglass 

En 

141. .34  GPa  (20.,50  Msi) 

21.99  GPa  (3.19  Msi) 

E-. 

=  £,< 

11.51  GPa  (1.07  Msi) 

21.99  GPa  (3,19  Msi) 

Go. 

=  Go 

0.(K)  GPa  (0.87  Msi) 

3.93  GPa  (0..57  Msi) 

G-. 

3  45  GPa  (0  ,50  Msi) 

2.90  GPa  (0.4.3  Msi) 

Vi: 

=  v,» 

0.20 

0  II 

Vs, 

0.43 

0,42 

TABLE  3  — 

Failure  loads  of  quasi-isotropic  short-heam  specimens. 

Test 

Load 

Shear  Strength 

1 

4.5.59  N  ( 1025  lbs) 

47091  kPa  (08,30  psi) 

2 

4982  N  (1120  lbs) 

5.5820  kPa  (8090  psi) 

3 

4093  N  ( 1055  lbs) 

46,^6<l  kPa  (6753  psi) 

4 

,5894  N  (1.325  lbs) 

5.301)0  kPa  (77(X,  psi) 

the  crack  front  is  sufficient  to  overcome  the  material  s  fracture 
resistance  or  toughness. 

Finite-element  method  (FE.V1)  anahsis  of  the  >pecimen  was 
performed  using  two  codes:  Engineering  Analysis  Language 
(EAL),  a  linear  code,  and  Cjeometrically  and  Materially  Non¬ 
linear  Analysis  of  Structures  (GAMNAS).  a  NAS.A-developed 
code  described  in  Ref  4.  Both  codes  were  run  using  meshes 
consisting  entirely  of  four-noded.  rectangular  elements,  using  a 
two-dimensional  (2-D)  plane-stress  model.  Both  meshes  consid¬ 
ered  one  ply  to  be  one  element  thick.  The  I:AL  mesh  used 
elements  of  aspect  ratio  5:1.  except  within  five-ply  thicknesses 
of  the  load  nose  and  support,  where  the  ratio  was  1:1.  These 
ratios  were  not  the  same  in  the  fiberglass  plies,  which  were  also 
''ne  element  thick  but  had  a  greater  real  thickness  than  the  other 
plies.  The  EAL  mesh  modelled  only  one  half  of  the  specimen, 
taking  advantage  of  the  symmetry  of  the  problem  about  the 
loading  axis.  Figure  2  shows  the  mesh  used  for  the  EAL  model. 
The  number  of  elements  is  1900,  and  the  total  number  of  degrees 
of  freedom  is  ,3940.  The  heavy  horizontal  line  in  the  figure  cor¬ 
responds  to  an  assumed  crack. 

Subsequent  EAL  models  involving  assumed  cracks  in  the  spec¬ 
imen  near  the  free  end  also  had  mesh  refinement  to  an  aspect 
ratio  of  1 : 1  in  the  vicinity  of  the  crack  tip.  These  models  were 
used  to  find  the  strain  energy  release  rates  associated  with  the 
assumed  cracks  using  the  zero-length  element  method  outlined 
in  Ref 


The  GAMNAS  model  did  not  use  a  symmetry  plane,  and  all 
elements  were  of  approximate  aspect  ratio  4:1,  except  for  ele¬ 
ments  in  the  fiberglass  plies,  which  were  of  approximate  aspect 
ratio  1,2:1.  There  are  2600  elements  and  5404  degrees  of  free¬ 
dom  in  this  mesh.  The  mesh  is  coarse  in  comparison  with  the 
EAL  discretization,  since  the  element  is  based  upon  a  nonlinear 
formulation. 

The  GAMNAS  model  was  used  to  predict  the  stress  field  in 
the  specimen  using  an  assumption  of  geometric  nonlinearity.  The 
EAL  model  was  used  to  predict  the  strain  energy  release  rate 
using  the  zero-length  element,  as  well  as  the  stress  field.  Com¬ 
parison  of  the  stress  field  results  from  the  two  codes  showed 
geometric  nonlinearity  effects  to  be  negligible.  This  supports  the 
use  of  the  linear  EAL  model,  which  was  much  simpler  to  use 
for  a  large  number  of  similar  but  nonidentical  configurations. 

Results  of  Numerical  Analysis 

Stress  Maps — Table  4  shows  the  locations  of  delamination  fail¬ 
ure  observed  visually  by  the  naked  eye  and  using  an  optical 
microscope  at  x  30  magnification  for  a  representative  group  of 
samples.  Figure  3  shows  the  ply-  ana  interface-numbering  scheme. 
The  table  lists,  for  each  interface,  the  number  of  specimens  in 


TABLE  4 — Locations  of  failure  of 
short-beam  specimens. 


Location, 

Number  of 

Interface  No. 

Failures 

14 

s 

11 

2 

7 

6 

1 

12 

1 

FIG.  3 — Ply  and  interface  niimherinft 
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which  final  failure  was  observed  for  that  interface.  Since  the 
layup  includes  50  plies,  there  are  49  interfaces,  of  which  Interface 
25  is  the  midplane.  It  is  seen  that  the  failures  took  place  well 
away  from  the  midplane.  Note  the  large  amount  of  variation  in 
the  locations  of  delamination.  No  single  interface  accounts  for 
even  a  majority  of  the  observed  failures.  This  would  not  be  the 
case  if  fracture  was  driven  by  a  simple  mechanical  process.  In 
other  materials,  such  variation  in  fracture  behavior  might  be 
attributed  to  the  random  distribution  of  material  defects.  Several 
types  of  processing  defects,  such  as  voids,  resin-rich  areas,  and 
bridged  fibers,  can  affect  fracture  behavior.  However,  such  de¬ 
fects  usually  do  not  drive  fracture  behavior  to  the  extent  found 
in  Table  4.  Also  note  that  several  specimens  showed  additional, 
secondary  cracking  above  the  midplane,  indicating  a  complex 
failure  process.  It  was  observed  that  the  delaminations  did  not 
extend  for  the  full  length  of  the  specimen.  Finally,  note  that  all 
of  the  specimens  in  which  one  primary  delamination  appeared 
to  the  naked  eye  had  that  single  delamination  below  the  mid¬ 
plane.  Figure  4  shows  actual  failed  specimens,  illustrating  the 
variability  in  crack  location.  The  locations  of  the  roller  supports 
and  applied  load  are  shown  schematically  in  Fig.  5. 


FIG.  4 — Failed  short  beam  specimens. 


FIG.  ?  —  Plotted  rcyion  for  stress  and  strain  maps. 
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The  GAMN  AS  analysis  yielded  stress  prrrfiles  for  the  speeimen 
under  load,  which  are  shown  as  Figs.  6,  7,  and  8  On  the  basis 
of  the  data  shown  in  Table  3,  a  load  of  4448  N  (UKK)  lbs)  was 
selected  for  these  calculations.  These  figures  show  one  half  of 
the  y-c  plane.  The  coordinate  system  can  be  understood  in  terms 
of  Fig.  1,  which  also  shows  the  y-z  plane.  The  origin  of  the 
coordinate  system  is  at  the  midplane  directly  beneath  the  load 
nose.  Thus,  the  load  is  applied  at  the  point  (0.  /2),  and  the 

supports  are  located  at  the  points  (±c,  A  schematic 

representation  of  this  region  is  shown  in  Fig.  5. 

Figure  6  maps  the  peel  stress  a  ,.  The  figure  shows  high  values 
of  compressive  stress  near  the  load  nose  and  the  support  and 
relatively  small  values  of  peel  stress  elsewhere. 


FIG.  b — Mti/>  aj  pt’t’l  slrcs.'i. 


Figure  7  maps  the  chordwisc  stress  ir  This  figure  is  much  as 
would  be  expected,  with  a  tspical  bending-stress  profile  and  the 
largest  magnitudes  near  the  load  nose.  As  a  check  for  reason¬ 
ableness,  a  two-dimensional  plot  of  ir,,  appears  on  the  insert  to 
Fig.  7.  These  values  are  taken  at  ^  7.9,37.S  m;n  |l)..M2.s  in.), 
a  location  well  removed  from  the  stress  concentrations  .it  the 
load  nose  and  support,  and  are  shown  as  the  light  line  on  the 
plot.  The  dashed  line  graphs  the  (r„  values  predicted  from  iso¬ 
tropic  beam  theory  for  a  homogeneous  be;im  with  smeared  prop¬ 
erties.  This  simple  model  well  represent;  the  m;ignitude  and 
trend  of  the  finite  element  data,  as  it  should. 

The  map  of  the  shear  stress  a,  .  shown  in  Fig.  8,  is  the  most 
significant.  Major  deviations  from  the  parabolic  stress  profile 
expected  from  beam  theory  are  seen.  High  values  of  shear  stress 
are  found  near  the  load  nose  and  support.  The  maximum  stress 
values  between  the  load  nose  and  support  arc  consistent  with 
the  predictions  of  beam  theory.  This  is  shown  by  an  insert  plot 
similar  to  that  of  Fig.  7. 

In  summary,  the  beam-theory  analysis  used  to  derive  Fq  I 
significantly  underestimates  the  magnitude  of  the  stress  beneath 
the  load  nose  and  near  the  support,  while  adequately  predicting 
it  elsewhere. 

Failure  /Vedir  tion.v  — Figures  6.  7,  and  8  indicate  that  it  is  un¬ 
likely  that  damage  initiates  at  the  specimen  free  ends,  which  are 
nearly  stress  free.  To  confirm  this,  failure  of  the  specimen  was 
analyzed,  using  the  maximum-strain  criterion.  Stress  data  from 
each  element  in  the  GAMNAS  model  were  converted  to  strain 
data,  and  these  strains  were  compared  to  the  materiiil  strain  to 
failure.  The  strains  to  failure  for  the  fiber,  matrix,  and  shear 
orientations  are  11  100.  6800,  and  20  000  gs,  respectively,  lire 
fraction  of  the  critical  strain  in  the  fiber,  matrix,  and  shear  ori¬ 
entations  was  ciilculated.  and  the  largest  of  these  was  selected, 
with  the  assigned  data  values  capped  at  100' f  of  strtiin  to  failure. 
The  resulting  data  provide  a  map  of  nearness  to  f..ilure  at  each 
location  in  the  specimen  and  are  shown  in  Fig.  0.  The  load  used 
for  this  figure  is  .S.338  N  ( 1200  lbs).  This  is  similar  to  the  actual 
failure  loads. 


FIG  7  —  Map  oj  chordwise  stress. 


FIG-  8  —  Map  oj  shear  stress. 


FIG.  9 — Map  of  muximuni  Inn  lioii  of  criliail  •noon 


FIG.  H) — Map  of  fraclion  of  cnlical  Mrain  as  a  rcsall  ol  chortini.sc 
stresses. 


Figures  1(>.  11.  and  12  show  nearness  to  failure  maps  for  the 
three  stress  components:  chordwise.  peel,  and  shear.  Considered 
along,  with  Fig.  9.  they  offer  insights  into  the  possible  causes  of 
failure  at  various  locations.  For  example.  Fig.  9  shows  that  failure 
is  imminent  in  the  region  near  the  load  nose;  reference  to  Figs. 
10,  11,  and  12  shows  that  this  is  primarily  due  to  chordwise 
stresses,  although  peel  and  shear  stresses  could  also  contribute. 
Failure  as  a  result  of  crushing  is  also  incipient  near  the  support, 
and  failure  as  a  result  of  chordwise  stresses  below  the  midplane 
near  the  load  nose  is  also  a  possibility.  In  contrast,  the  midplane, 
where  beam  theory  predicts  the  largest  shear  stresses,  is  the 
region  least  vulnerable  to  failure.  This  is  consistent  with  the 
observed  cracks  in  the  tested  specimens,  which  failed  below  the 
midplane  and  sometimes  above  H,  b'.it  never  on  or  near  the 
midplane. 


FORTSON  AND  ARMANIOS  ON  DAMAGE  PROCESSES  253 


FIG.  1 1  —  Map  of  fraction  of  critical  strain  as  a  resiili  ol  peel  stress 


Flfj  12  —  Map  of  fraction  of  critical  strain  as  a  result  of  shear  stresses. 


The  maximum  strain-to-failure  map  of  Fig.  9  was  modified  to 
include  additional  failure  modes  observed  by  McCleskey'  for 
laminates  made  of  glass-epoxy  material.  McCIcskey's  laminates 
consisted  of  eight  plies,  each  of  thickness  ().82.‘>.s  mm  (()0.12.s 
in.),  and  were  of  layup  (9U/-45'45.'0l,.  When  these  specimens 
were  tested  in  three-point  bending,  a  ()“  pl\  near  the  midplane 
was  observed  to  crack  at  an  angle  of  4.“''  to  the  midpuine,  as 
shown  in  Fig.  1.4,  McCleskey  suggested  that  these  failures  were 
due  to  the  state  of  nearly  pure  shear  near  the  specimen  midplanc 
resolving  itself  into  a  tensile  stress  at  a  45°  angle  to  the  midplane. 
This  stress  would  act  upon  the  matrix  of  a  tl°  ply. 

'S.  F.  McCleskcv,  private  communication.  .^Ist  A1/\A .ASMF-  .AS(  H 
AHS.  ASC  Structures.  Structural  Dynamics  and  Mater'als  Conlerence. 
Long  Beach.  CA.  2-4  April  149(1, 
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This  effect  could  be  modelled  as  shown  in  Fig.  14,  It  is  assumed 
that  any  normal  stress  in  the  y-z  plane  leads  to  matrix-driven 
failure  in  a  0°  ply.  Thus,  to  determine  the  maximum  effect  of 
normal  stres.ses  upon  such  a  ply,  it  is  necessary  to  determine  the 
largest  normal  stress  operating  in  the  y-z  plane;  if  a,,  is  assumed 
to  be  zero,  then  this  largest  normal  stress  is  the  principal  stress 
with  the  largest  magnitude.  The  percentage  of  strain  to  failure 
can  then  be  calculated  using  this  largest  principal  stress,  and  the 
appropriate  tensile  or  compressive  modulus  for  loads  applied 
normal  to  the  fibers. 

A  similar  phenomenon  might  take  place  in  the  45"  plies,  lead¬ 
ing  to  a  failure  surface  as  shown  in  Fig.  15.  To  predict  this  second 
additional  failure  mode,  it  would  be  necessary  to  determine  the 
traction  normal  to  the  fracture  surface  shown  in  this  figure,  and 
therefore,  to  perform  a  three-dimensional  stress  transformation, 
assuming  the  out-of-plane  stresses  to  be  zero.  The  percentage  of 
strain  to  failure  for  this  fracture  mode  could  then  be  calculated 
using  the  traction  resolved  on  this  face,  and  the  appropriate 
tensile  or  compressive  modulus  for  loads  applied  normal  to  the 
fibers. 

The  analysis  that  produced  Fig.  9  was  modified  to  include  the 
two  additional  failure  modes  and  was  repeated,  with  results  as 
shown  in  Figs.  16,  17.  and  18.  The  areas  in  which  Figs.  9  and  16 
differ  are,  therefore,  areas  in  which  one  of  the  McCleskey  failure 
modes  determines  fracture.  Comparison  of  Figs.  9  and  16  shows 
one  change  as  a  result  of  these  new  modes  of  failure;  some 
interior  plies  become  much  more  highly  strained  in  the  region 
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ZERO- DEGREE  FIBER 


PRINCIPAL  STRESS 

FIG  14 — Moilelliny  of  first  McC  h'ski's  tiiilurc  modi' 


between  the  load  nose  and  the  support.  These  are  U  plies,  with 
the  new  failure  mode  being  the  one  shown  in  Fig  l.s.  Figures 
1 7  and  18  show  the  contributions  of  the  two  individual  Mcf  leskev 
failure  modes.  In  summary.  Fig.  17  includes  the  two  failure  modes 
observed  by  McC'leskey,  By  comnarison  with  Fig.  9,  the  map  of 
Fig.  16  shows  a  significantly  increased  ne.irness  to  failure  iii  the 
interior  of  this  specimen.  While  there  is  not  .  nough  of  an  increase 
to  make  this  a  principal  initial  failure  location,  the  possibililv  of 
cracking  in  the  interior  of  the  specimen  is  significantly  enhanced. 

Strain  Etwry,y  Relea.se  Rale  AnuA  st' —  f;AI.  modelling  was  used 
to  construct  strain  energy  release  rate  profiles  at  three  locations 
In  the  first  profile,  a  crack  of  length  1 .778  mm  |(MI7  in  ).  centered 
beneath  the  load  nose,  was  considered.  It  was  necessarv  to  make 
the  crack  this  length  to  distance  the  crack  tip  from  the  region  of 
the  symmetry  plane  used  in  the  model.  A  schematic  of  the  mesh 
is  shown  in  Fig,  2  for  a  crack  between  the  45  -45  plies  at  In¬ 
terface  .52.  This  crack  was  assumed  at  each  of  the  various  inter¬ 
faces.  and  the  resulting  strain  energy  release  rates  were  calcu- 

iZ- 


45"  FRACTURE 


FKj.  15 — The  wcomi  McC  leskcv  I'.uliin"  nioilc. 


FIG.  — \fup  f>/  nvarncw  u>  lailiirv  int  <>i poriinni^  \tf(  lt  \kc\  Uiilure 

nuides. 
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FIG  17  —  Map  i>l  framoii  of  cnlical  strain  as  a  result  ul  first  Me(  leskev 
failure  mode 


FIG.  IK — Map  Ilf  Inieiioii  of  t  ritual  strain  as  a  result  ot  seeond  Me(  leskes 
fatltire  mode. 

lated.  This  addressed  specimen  behavior  in  the  region  of  maximum 
absolute  chord  stress  and  peel  stress.  In  the  second  profile,  this 
was  repeated  for  an  assumed  crack  of  length  (I,2,s4  mm  (fl.Ol 
in.),  centered  10.16  mm  (0.40  in.)  from  the  load  nose,  near  the 
region  of  maximum  absolute  shear  stress.  Finally,  a  crack  of 
length  2..S4  mm  (0. 10  in. )  was  assumed  at  the  free  end.  suggested 
by  initial  observation  to  be  the  site  of  crack  initiation.  The  applied 
load  in  all  cases  was  4448  N  ( KKK)  lbs),  as  with  the  stress  maps. 

The  profile  for  the  assumed  crack  near  the  free  end  showed 
what  Figs.  6,  7,  and  8  suggest:  negligible  values  of  the  strain 
energy  release  rate,  too  small  to  produce  failure.  The  data  for 
the  first  two  profiles  appear  as  Fig.  I V.  The  distance  y  is  as  defined 
in  Figs.  6,  7,  and  8. 

Note  that  all  data  points  shown  in  Fig.  14  are  taken  from 
locations  at  least  1.0 '6  mm  (0.04  in.)  from  the  surface  of  the 
specimen.  Thus,  for  ail  data  points,  the  peel  stress  is  negligible, 
as  indicated  by  Fig.  6.  Consequently,  crack  growth  is  predomi¬ 
nantly  Mode  11.  The  spikes  in  the  Ci  distribution  result  from 
variations  in  elastic  properties  from  ply  to  ply. 


FIG  ly  —  Strum  eneri;\  release  rate  data  trom  I.AI. 

As  shown  in  Fig.  14.  the  strain  energy  release  rate  for  assumed 
interlaminar  cracks  beneath  the  load  nose  increases  rapidly  as 
the  flaw  is  placed  nearer  to  the  point  of  application  of  the  load. 
However,  for  flaws  assumed  in  the  range  of  high  shear  stress 
away  from  the  load  nose,  strain  energy  release  rates  remain 
relatively  negligible.  Figure  14  shows  that  cracking  can  occur  in 
plies  near  the  load  nose,  but  is  unlikely  away  from  this  region. 
This  supports  the  data  of  Fig.  4.  which  indicate  that  overstress 
failure  is  found  in  this  same  region. 

Results  of  Fraaographic  Analysis 

Microscopic  analysis  of  a  failed  specimen  showed  several  per¬ 
tinent  features,  many  of  which  were  visible  with  the  naked  eye 
after  the  specimen  had  been  sputter-coated.  Directly  beneath 
the  load  nose,  many  cracks  were  observed  running  perpendicular 
to  the  direction  of  load  application.  One  such  is  shown  in  Fig. 
21).  This  crack  was  approximately  b  .I.K  mm  (0  2.“'  in.)  long  and 
was  found  beneath  the  load  nose,  above  the  midplane. 

Figure  21  shows  some  smaller  flaws  in  the  same  1  he 

presence  of  such  flaws  supports  the  analytical  determination  that 
extensive  damage  as  a  result  of  crushing  in  this  region  will  lead 
to  the  growth  of  cracks. 

Figure  22  shows  the  failure  of  a  4.s-  ply  below  the  midpl.me. 
near  the  load  nose.  This  ply  has  experieneed  matrix  failure, 
splitting  in  two  along  the  fiber  direction.  This  supports  the  an¬ 
alytical  prediction  of  ply  failure  in  this  location  as  a  result  of 
chordwise  stresses. 

The  presence  of  matrix  failure  suggests  an  explanation  for  the 
failure  of  all  specimens  below  the  midplane.  In  this  region,  the 
chordwise  stresses  are  positive,  leading  to  possible  tensile  failure 
of  the  matrix,  which  is  more  likely  to  lead  to  eompiete  speeimen 
failure  than  is  the  compressixe  failure  which  is  to  be  found  aboxe 
the  midplane. 

Conclusion 

Finite-element  analyses  have  been  used  to  describe  the  be¬ 
havior  of  a  quasi-isotropic  short  beam  specimen.  Short  beam 
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FICi  2n  (  raik  hi’m-aih  load  n<>\f 


in. 


FTCi  21  —  SituiH  i  riiik\  hi'iii'iilh  Imid  >iii\i' 


specimens  hiive  hccn  tested  and  the  results  eumpareil  with  the 
analyses,  and  the  fracture  surfaces  of  failed  specimens  base  been 
examined. 

Fractographic  analvses  show  th.it  numerous  sm.ill  cracks  arc 
created  directly  beneath  the  load  nose.  possibK  c.iusing  fracture 
to  proceed  from  this  location  This  is  supported  b\  the  presence 
of  large  magnitudes  of  peel,  shear,  and  chordwise  stresses  m  the 
vicinity  of  the  load  nose,  as  shown  by  the  1  l.M  an.ilssis.  and  b\ 
comparison  of  strain  energy  release  rate  values  obtained  at  the 
same  locations. 

The  following  mechanism  may  be  postulated.  Stresses  c.iuse 
matrix  failures  beneath  the  load  nose  Above  the  midplane,  this 
failure  is  caused  by  crushing  stresses  in  the  c  direction;  below 


FKi.  22-  \liliri\  laihirc  nl  n  4'  pl\ 


the  midpliine.  it  is  c.iused  by  chordwise  tensile  stresses  These 
regions  of  damage  are  sites  from  which  dclamin.ition  cr.icks  c.in 
propagate.  Fventu.illy  .  the  region  benc.iih  the  midpl.me  becomes 
sufficiently  damaged  to  produce  a  del.imination  th.it  will  prop¬ 
agate  to  failure.  The  cr.icks  origiiiiiting  .ibove  the  midplane  .ire 
subjected  to  a  compressive  clurrdwisc  stress  and  to  a  larger  com¬ 
pressive  peel  Stress. 

The  strength  of  materials  method  neglects  the  chordwise  and 
crushing  loads  which  are  the  driving  forces  causing  f.iilure  in  this 
type  of  specimen.  Therefore,  the  three-point  short-beam  shear 
test  does  not  measure  the  interlaminar  shear  strength  of  a  quasi- 
isotropic  l.iminated  composite  Because  of  the  inter.iction  of 
damage  modes  in  this  configuration,  the  vjuantitativ e  vlata  c.innot 
be  interpreted  as  a  me.isure  of  .1  simple,  isolatcrl  tailure  mode. 
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